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Procede de determination de parametres de fonctionnement d'une 
installation comprenant une enceinte de refroidissement. 

La presente invention est relative a un procede de determination de 
parametres de fonctionnement d'une installation de refroidissement thermique 
d'articles. 

5 ^invention s'applique notamment aux installations de surgelation 

d'articles alimentaires. 

Les installations connues de surgelation comportent, par exemple, une 
enceinte ou tunnel, de surgelation traverse de part en part par un convoyeur a 
bandes sur lequel sont deposes des articles a congeler, le convoyeur circulant en 
10 continu ou sequentiellement au travers du tunnel de surgelation. 

Une enceinte cryogenique utilise un fluide inerte a basse temperature 
qui echange de la chaleur directement par contact avec les articles a surgeler. 

De maniere classique, une enceinte cryogenique utilise soit la neige 
carbonique (-80°C), I'air liquide ou I'azote liquide (-196°C) comme vecteur de 
15 froid. La neige carbonique permet le transport de produits frais ou surgeies sans 
crainte de rupture de la chatne de froid. L'azote et Fair liquide permettent soit la 
surgelation individualisee de produits alimentaires, soit le durcissement de pro- 
duits fragiles, deformables ou collants (du type creme glacee ...). 

Les parametres de fonctionnement des installations forment des recet- 
20 tes creees de maniere experimental. Une recette memorise les parametres de 
reglage d'une installation pour une production donnee. 

II n'existe pas aujourd'hui de capteur capable de mesurer en continu et 
sans contact, la temperature interne d'un article, de sorte que la determination 
des recettes passe par des procedes de determination de parametres de fonc- 
25 tionnement qui sont ensuite ajustes a I'aide de tests realises en situation, les- 
quels tests sont le plus souvent destructifs. 

Les procedes de determination de recettes comprennent generale- 
ment une etape de determination d f une consigne portant sur la temperature de 
sortie des articles suivie d'une etape de determination de parametres de fonc- 
30 tionnement initiaux puis la reiteration d'un cycle de test comprenant une etape de 
prediction de la temperature de sortie des articles, une etape de comparaison de 
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la temperature predite avec la consigne et, en cas de difference, une etape de 
modification des parametres de fonctionnement. 

Le cycle de test est renouvele jusqu'a ce que les parametres de fonc- 
tionnement permettent d'obtenir une temperature predite sensiblement proche de 
5 la temperature de consigne. 

II est a noter que toutes ces etapes sont deduites et realisees par 
I'operateur de fagon experimental, en tenant compte de son savoir faire, son 
experience, de ses connaissances. 

Si Ton examine le systeme constitue de I'enceinte et de la charge d'ar- 
10 tides, plusieurs parametres peuvent influer sur la temperature de Particle en sor- 
tie : le debit de production qui, pour un taux de chargement donne, implique une 
variation du temps de sejour dans I'enceinte, le debit du fluide qui agit sur le profil 
des temperatures, la temperature d'entree des articles, le profil convectif de 
I'enceinte, et le taux de chargement. 
15 Le systeme est done un systeme multivariates et les procedes de de- 

termination de parametres existants, et notamment Tetape de prediction de la 
temperature des articles en sortie, ne prennent pas en compte ces elements. 

Dans les etapes de prediction de I'etat de la technique, afin de traiter 
des systemes monovariables, on a ete amene a considerer le profil convectif et le 
20 taux de chargement comme constants, et a fixer le debit de production, la tempe- 
rature des articles en entree de I'enceinte de surgelation ainsi que d'autres 
parametres de fonctionnement des installations. 

De ce fait, les recettes determinees par les procedes de determination 
existants sont relativement imprecises et necessitent la mise en oeuvre d*une 
25 production suivie de tests destructifs sur les articles. 

La presente invention vise a remedier a ce probleme en definissant un 
procede de determination de parametres de fonctionnement, precis et aisement 
mis en oeuvre. 

A cet effet, ['invention a pour objet un procede de determination des 
30 parametres de fonctionnement d'une installation de refroidissement thermique 
d'articles, comprenant une enceinte traversee par lesdits articles d f une entree a 
une sortie et utilisant un fluide de refroidissement, le procede comprenant : 

- une etape de determination d'une consigne de temperature d'articles 
en sortie de ladite enceinte ; 
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- ladite etape de prediction du comportement desdits articles est reali- 
see en outre, a partir de parametres de fonctionnement de ladite installation ; 

- lesdits parametres de fonctionnement de ladite installation compren- 
nent la temperature desdits articles en entree de ladite enceinte ; 

5 - ladite etape de prediction du comportement desdits articles est opti- 

misee par des calculs de modification dudit maillage desdits articles selon des 
suites mathematiques ; 

- ladite etape de prediction du comportement desdits articles est opti- 
misee par suppression des calculs de prediction pour des points spatiaux et tem- 

10 porels dudit maillage desdits articles pour lesquels les variations d'enthalpie sont 
inferieures a un seuil predetermine ; 

- ladite etape de prediction de la temperature desdits articles en sortie 
de ladite enceinte est fondee sur ladite etape de prediction du comportement de 
ladite enceinte ainsi que sur ladite etape de prediction du comportement desdits 

15 articles; 

- ladite etape de modification des parametres de fonctionnement com- 
prend une etape de modification manuelle d'au moins une partie des parametres 
de fonctionnement ; 

- ladite etape de modification des parametres de fonctionnement com- 
20 prend la modification automatique d'au moins une partie desdits parametres de 

fonctionnement ; et 

- ladite etape de modification des parametres de fonctionnement com- 
prend la modification d'au moins un des parametres choisis par le groupe consis- 
tant en : 

25 - le debit dudit fluide de refroidissement ; 

- le temps de sejour desdits articles dans ladite enceinte ; 

- le debit de gaz extrait de ladite enceinte ; 

- la mise en vitesse des gaz ; 

- la recirculation des gaz ; et 

30 - la balance entre les entrees d'air et les sorties de gaz. 

[..'invention sera mieux comprise a la lecture de la description qui va 
suivre, donnee uniquement a titre d'exemple et faite en se referant aux dessins 
annexes, sur lesquels : 
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- une etape de determination de parametres de fonctionnement initiaux 
pourladite installation ; et 

- un cycle de test des parametres de fonctionnement comprenant : 

- une etape de prediction de la temperature d'articles en sortie de 

5 ladite enceinte ; 

- une etape de comparaison de la temperature de consigne avec 
la temperature predite ; et 

- si ladite etape de comparaison revele un ecart superieur a un 
seuil predetermine, une etape de modifications de parametres de fonctionnement 

10 de ladite installation et une reiteration du cycle de test, 

ladite etape de prediction etant realisee a partir de parametres de 
fonctionnement de ladite enceinte, de caracteristiques thermodynamiques et 
physiques de ladite enceinte et de caracteristiques thermodynamiques et physi- 
ques desdits articles. 

1 5 Suivant d'autres caracteristiques : 

- ladite etape de prediction comporte une etape de prediction du com- 
portement de ladite enceinte fondee sur la resolution de bilans thermiques sur 
des tranches elementaires du volume de ladite enceinte, realisee au moins a par- 
tir de caracteristiques thermodynamiques dudit fluide de refroidissement et de 

20 caracteristiques thermodynamiques et physiques de ladite enceinte ; 

- ladite etape de prediction du comportement de ladite enceinte est re- 
alisee en outre, a partir de parametres de fonctionnement de ladite installation ; 

- lesdits parametres de fonctionnement de ladite installation represen- 
ted au moins Tun des elements choisis dans le groupe consistant en : 

25 - la vitesse d'un convoyeur de transport desdits articles au tra- 

vers de ladite enceinte ; 

- le taux de chargement ; et 

- la ventilation de Tatmosphere de ladite enceinte ; 

- ladite etape de prediction comporte une etape de prediction du com- 
30 portement desdits articles fondee sur la resolution de Tequation de conservation 

de la chaleur discretisee et appliquee a un reseau de points spatiaux et tempo- 
rels constituant un maillage desdits articles, realisee au moins a partir de caracte- 
ristiques thermodynamiques et physiques desdits articles ; 
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- la Fig.1 represente un schema synoptique illustrant une installation 
de refroidissement ; 

- la Fig.2 est un organigramme general du procede de I'invention ; 

- la Fig.3 illustre la modelisation numerique des articles a traiter ; 

5 - la Fig.4 illustre la modelisation numerique de I'enceinte de refroidis- 

sement ; et 

- la Fig. 5 represente I'organigramme detaille du cycle de test du pro- 
cede de Tinvention. 

Sur la figure 1, on a represente une installation classique de traitement 
10 d'articles alimentaires, pour laquelle des parametres de fonctionnement sont de- 
termines par un procede selon Tinvention. 

Cette installation comprend une enceinte ou tunnel cryogenique 2, de 
type classique, permettant la congelation d'articles alimentaires P par leur mise 
en presence avec un fluide cryogenique 4 vehicule par une ligne d'alimentation 5, 
15 a partir d'une source quelconque. 

Par exemple, I'enceinte 2 a une forme de parallelepipede rectangle. 
Comme cela est dit precedemment, le fluide cryogenique 4 utilise peut 
etre par exemple, de la neige carbonique ou de I'azote liquide et est injecte en un 
ou plusieurs endroits de I'enceinte 2. 
20 Cette enceinte 2 est associee a un convoyeur 6 de type classique, 

permettant ('introduction des articles P dans I'enceinte 2 et leur extraction et fonc- 
tionnant soit de maniere sequentielle, soit de maniere continue. 

L'installation presente plusieurs parametres de fonctionnement, a sa- 
voir le profil des temperatures dans I'enceinte, le temps de sejour des articles P 
25 dans I'enceinte 2 ou la vitesse de deroulement du convoyeur 6, et la temperature 
d'entree des articles P. 

L'installation comporte enfin des moyens de commande 12 de la quan- 
tity de fluide cryogenjque 4 injectee dans I'enceinte 2. 

Ces moyens 12 comportent des moyens de pilotage 14 du debit de 
30 fluide cryogenique 4. Par exemple, les moyens de pilotage 14 sont constitues par 
des systemes d'electrovannes ou des vannes proportionnelles de type classique, 
disposes sur la ligne d'alimentation 5 en fluide cryogenique 4. 
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Avantageusement, Installation comporte egalement un systeme de 
ventilation de gaz controlant les flux de gaz et la ventilation de I'atmosphere de 
I'enceinte 2. 

Par exemple, ce systeme est compose de ventilateurs specifiques 
5 permettant la mise en Vitesse des gaz, de ventilateurs controlant la recirculation 
des gaz et d'une combinaison entre des ventilateurs et des portes mobiles 
controlant la balance entre les entrees d'air et les sorties de gaz. 

En reference a la figure 2, on va maintenant decrire un organigramme 
general du procede de determination des parametres de fonctionnement selon 
10 I'invention. 

Ce procede debute par une etape 16 de saisie d'une consigne portant 
sur la temperature de sortie des articles apres refroidissement thermique. 

Uetape 16 est suivie d'une etape 18 de determination des parametres 
de fonctionnement initiaux. Les parametres determines au cours de cette etape 
15 18 sont des parametres connus tels que les caracteristiques mecaniques de Ten- 
ceinte 2 ou encore les caracteristiques physiques et thermodynamiques des arti- 
cles P et des parametres variables, tels que les parametres de fonctionnement 
de installation qui sont fixes arbitrairement 

Le procede comporte ensuite une etape 20 de prediction de la tempe- 
20 rature des articles P en sortie de I'enceinte 2. 

Cette etape 20 comporte une etape 22 de prediction du comportement 
de I'enceinte 2 et une etape 24 de prediction du comportement des articles P. 

Uetape 22 de prediction du comportement de Tenceinte 2, permet de 
predire par calcul ainsi que cela est decrit plus loin en reference a la figure 4, le 
25 profil theorique des temperatures du fluide cryogenique a I'interieur de I'enceinte 
2. 

Les resultats delivres par I'etape 22 dependent des caracteristiques 
thermodynamiques du fluide cryogenique 4, des caracteristiques convectives de 
I'enceinte 2 ainsi que les caracteristiques des moyens d'injection du fluide cryo- 
30 genique 4, des caracteristiques du systeme de ventilation et des caracteristiques 
physiques de Tenceinte 2. 

L'etape 22 prend egalement en compte des parametres de fonction- 
nement de 1'installation comme la vitesse du convoyeur 6. 
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L'etape 24 de prediction du comportent des articles P permet de de- 
terminer par calcul, ainsi que cela est decrit plus loin en reference a la figure 3, 
des variations d'enthalpie des articles P en fonction de leur environnement exte- 
rieur et de leur temperature initiale. 
5 Ainsi, les resultats delivres par l'etape 24 de prediction du comporte- 

ment des articles P, dependent de leurs caracteristiques physiques et thermody- 
namiques. 

Les etapes 22 de prediction du comportement de I'enceinte 24 et de 
prediction du comportement des articles P, sont couplees Tune a I'autre ainsi que 
10 cela est decrit plus en detail en reference a la figure 5, afin de delivrer une tem- 
perature theorique des articles P en sortie de I'enceinte 2. 

Ainsi, la prediction realisee au cours de I'etape 20 de prediction de la 
temperature des articles P en sortie de I'enceinte 2 prend en compte les caracte- 
ristiques thermodynamiques et physiques de I'enceinte 2 et des articles P, ainsi 
15 que les parametres de fonctionnement de Installation. 

De ce fait, la determination de la temperature des articles P en sortie 
de I'enceinte 2 est dynamique, parametrable et d'une grande precision. 

L'etape 20 de prediction de la temperature des articles P en sortie de 
I'enceinte 2 est suivie d'une etape 26 de comparaison de la temperature predite 
20 lors de l'etape 20, avec la consigne de temperature determinee lors de I'etape 16. 

Par exemple, cette etape 26 de comparaison permet la prise en 
compte d'un intervalle de tolerance de I'ordre de quelques degres autour de la 
consigne de temperature de sortie determinee lors de l'etape 16. 

Quand I'ecart entre la temperature predite et la consigne est superieur 
25 a un seuil predetermine, l'etape 26 de comparaison est suivie d'une etape 28 de 
modification de parametres de fonctionnement de ('installation. 

Par exemple, les parametres de fonctionnement modifies au cours de 
cette etape 28, comprennent la vitesse du convoyeur 6 et les parametres interve- 
nant dans la determination du profil theorique des temperatures du fluide 4 dans 
30 I'enceinte 2, soit par exemple le debit du fluide cryogenique 4, le controle de la 
ventilation a I'interieur de I'enceinte 2 et le taux de chargement du convoyeur 6. 

Les modifications des parametres de fonctionnement en cours de 
l'etape 28 peuvent etre realisees directement par un operateur ou etre faites de 



WO 2004/092667 



PCT/FR2004/050121 



8 

maniere automatique par un calculateur prenant en compte des limites maximale 
et minimale pour chacun des parametres. 

De plus, les modifications peuvent affecter un ou plusieurs parametre 
a chaque fois. 

5 II est egalement possible de definir un ordre de modification des para- 

metres de fonctionnement pour tenter d'atteindre la temperature de consigne en 
ne modifiant qu'un parametre a la fois. Si la consigne ne peut etre atteinte en 
modifiant un premier parametre entre des valeurs limites. Ce parametre est fixe 
sur une valeur limite ou une valeur moyenne et aux iterations suivantes, le para- 

10 metre suivant dans la liste est modifie. 

A Tissue de Tetape 28 de modification des parametres de fonctionne- 
ment, le procede reprend a Tetape 20 de prediction de la temperature des articles 
P en sortie de Tenceinte 2, le procede formant ainsi un cycle de test des parame- 
tres de fonctionnement de Installation, comprenant Tetape 20 de prediction de la 

15 temperature des articles en sortie, Tetape 26 de comparaison de la temperature 
predite avec la temperature de consigne et Tetape 28 de modification des para- 
metres de fonctionnement. 

Ce cycle de test, designe par la reference numerique generate 29, est 
repete jusqu'a ce que la comparaison effectuee lors de Tetape 26 entre la tempe- 

20 rature predite et la temperature de consigne revele a un ecart inferieur a une va- 
leur de seuil predeterminee. 

Le cycle 29 est alors interrompu et Tetape 26 est suivie d'une etape 30 
d'enregistrement des derniers parametres de fonctionnement testes, qui forment 
ainsi une recette. 

25 La prediction realisee au cours de Tetape 20 etant d'une grande preci- 

sion et prenant en compte des parametres de fonctionnement de Tenceinte 2, 
des caracteristiques thermodynamiques et physiques de Tenceinte 2 et des ca- 
racteristiques thermodynamiques et physiques de Tarticle P, la recette determi- 
nee par le procede de Tinvention est precise et proche du fonctionnement reel. 

30 De plus, une telle recette est facilement adaptee a des modifications 

des conditions de fonctionnement. 

Par exemple, si un parametre de fonctionnement tel que la tempera- 
ture d'entree des articles dans Tenceinte 2 est modifie, la recette peut etre corrig- 
ee par ('execution du procede de Tinvention en utilisant les parametres actuels de 
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recette comme parametres de fonctionnement initiaux lors de I'etape 18, I'execu- 
tion du procede de invention permettant de determiner rapidement et simple- 
ment les corrections a apporter aux parametres de fonctionnement pour obtenir 
la bonne temperature de sortie des articles P. 
5 En reference aux figures 3, 4 et 5, on va maintenant decrire plus en 

detail I'etape 20 de prediction de la temperature des articles P en sortie de Ten- 
ceinte 2. 

Sur la figure 3, on a represents un exemple de maillage d'un article 
alimentaire P tel que mis en ceuvre au cours de I'etape 24 de prediction du com- 
1 0 portement des articles P. 

La prise en compte, dans le procede de determination de parametres 
de fonctionnement, des caracteristiques thermodynamiques et physiques des 
articles P lors de I'etape 24 de prediction de leur comportement, est basee sur 
une modelisation des articles P a laquelle est appliquee I'equation de conserva- 
1 5 tion de la chaleur discretisee. 

En effet, I'equation de la conservation de la chaleur ne peut pas etre 
resolue en tout point de I'espace et en tout instant, par une fonction integrate 
simple. 

La methode employee consiste a discretiser cette equation de sorte 
20 que Ton ne la resout plus que sur des points spatiaux et temporels appeles 
noeuds et designes par la reference numerique generate 32. 

Apres definition d'un maillage de Tarticle P, on applique I'equation de 
conservation de la chaleur a chacun des nceuds 32. 

On obtient ainsi un systeme d'equation qu'il faut resoudre pour connaT- 
25 tre I'etat thermique de I'article P dans le temps et dans I'espace. 

X, Y et Z sont des axes definissant un repere spatial orthonorme au- 
tour de I'article P. T est la temperature de I'article P exprimee en kelvin (K), et C 
sa chaleur specifique exprimee en watt par kilogramme et par kelvin (W/(Kg*K)). 
30 Les articles alimentaires P que Ton surgele sont generalement consti- 

tues de corps differents. 
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Cela signifie que le changement de phase s'accompagne cTune varia- 
tion de temperature et que I'equation de la conservation de la chaleur peut tou- 
jours s'appliquer. 

En revanche, si Ton est amene a traiter un corps pur, I'equation n'est 
5 plus continue. Dans ce cas, on simplifie le probleme en modifiant la table 
d'enthalpie du corps pur pour que le changement de phase engendre une faible 
variation de temperature. 

La discretisation est faite grace a la methode mathematique des diffe- 
rences finies en regime variable. 
10 De maniere connue, celle-ci peut etre effectuee de deux manieres. 

La premiere, la discretisation implicite, a I'avantage d'etre stable quelle 
que soit la configuration spatiale et temporelle. A un instant donne, elle permet 
de determiner la temperature d'un noeud 32 en fonction de la temperature des 
noeuds voisins au meme instant. Elle implique cependant des conditions aux 
15 frontieres constantes et une resolution matricielle du systeme d'equation forme 
par chacun des noeuds 32. 

La seconde, la discretisation explicite, permet de determiner directe- 
ment le temperature d'un noeud 32 a un instant T+AT d'apres les conditions a 
llnstant T. Le resultat est immediat, en revanche, il faut choisir un pas de temps 
20 adapte de maniere a eviter Tinstabilite du models. 

La premiere methode est preconisee dans le cas ou Ton cherche a ob- 
tenir principalement la temperature de surface d'un article, ce qui correspond a 
I'operation communement appelee operation de « croutage ». La seconde est 
preconisee lorsque Ton veut faire de la surgelation et connaTtre la temperature a 
25 cceur d'un article. 

Le maillage de I'article P est un probleme crucial. De cela depend di- 
rectement la simplicite du traitement ulterieur et la precision des resultats. 

Un nombre de nceuds important amene une grande precision dans le 
resultat mais impose un temps de calcul eleve. II s'agit de trouver un compromis 
30 entre la precision et le temps de calcul. 

Par exemple, pour le cas d'un produit alimentaire de dimensions exte- 
rieures 100 x 60 x 10 mm avec un maillage regulier tous les millimetres, il faut 
plus de 17000 noeuds et autant d'equations pour definir le comportement de 
Tarticle P. 
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Dans la version decrite de invention, Tetape 24 comporte des calculs 
qui permettent d'optimiser le maillage qui est effectue pour la prediction du corrv 
portement des produits P. 

Par exemple, dans le cas du croutage, on surveille plus particuliere- 
5 ment la solidification d'une fine epaisseur de la peau du produit par changement 
de phase. II faut done un maillage dense en peripherie et plus large a coeur. 

Pour ne pas saisir manuellement les coordonnees de chacun des 
nceuds et pour conserver des relations simples entre les noeuds et faciliter le trai- 
tement, une solution consiste a distribuer des noeuds dans chaque direction de 
10 I'espace a I'aide par exemple d'une progression geometrique, comme cela est 
represents sur la figure 3. 

Par exemple, sur I'axe X, on distribue les noeuds de la maniere sui- 
vante : on considere Ax la valeur du premier terme qui correspond a I'abscisse du 
premier noeud, et r la raison, differente de 1, de la suite geometrique mis en oeu- 
15 vre. La valeur du n 6me terme est : Ax^r"" 1 , cela correspond a la position sur I'axe X 
du n 6me nceud. La somme des n premiers termes est : 

5 = Ax-f '— 

i-r 

La figure 3 represente le positionnement des noeuds selon ce maillage 
sur un article P parallelepipedique ou Ton a impose une condition de parite sur le 
20 nombre de noeuds de maniere a simplifier la resolution. 

On obtient ainsi une valeur correspondant a la dimension du produit P 
selon I'axe X : 



avec R = — — 
l + r 

Avec I qui correspond a I'abscisse du noeud central sur cette lon- 

25 gueur : 



Ln 

1= 



-XT 



A 1 R 

Ax 



L'imprecision sur I'axe X s'exprime alors de la maniere suivante : 
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Ln 

\I1<1 



Ax 



Ln 



<l/2 « ^-^ r '^-^^ 1 + lt r '))< ip+ ^ r '[ r i- 



Pour aboutir a des calculs simples on fixe les imprecisions sur les trois 
axes a une meme valeur. Cela induit une erreur sur les dimensions du produit P 
5 qui est acceptable dans le cas ou Ton s'interesse uniquement aux temperatures 
sur une faible epaisseur de peau et que les temperatures a coeur varient peu, 
comme cela est le cas dans les operations de croutage. 

Dans le cas des operations de surgelation ou Ton cherche a determi- 
ner la temperature a cceur du produit, un terme correctif peut etre insere dans les 
10 formules. Dans le cas de I'axe X, on insere le terme correctif suivant : 

L P 



Ar'= 



Une autre methode d'optimisation possible consiste a diminuer le 
temps de traitement en omettant certains calculs. 

En effet, sur chaque nosud, on additionne le flux thermique aux six fa- 
15 ces de son volume elementaire. Cependant, il existe des zones ou les effets 
thermiques sont assimilables a des problemes monodimensionnels. 

Pour exploiter cette particularite, on decompose le traitement en som- 
mant non plus globalement les flux thermiques sur chaque face, mais selon cha- 
que direction. Sur chaque direction, on resout, pour un pas de temps AT, les 
20 equations aux noeuds en allant de la frontiere vers le cceur, jusqu'a ce que la va- 
riation d'enthalpie soit consideree comme etant negligeabie car inferieure a un 
seuil predetermine. 

En effectuant cette operation dans chaque direction, on definit un vo- 
lume du produit P englobant tous les nceuds pour lesquels les variations 
25 d'enthalpie seront negligeables, et done pour lesquels aucun calcul ne sera fait. 

On peut ainsi economiser du temps de calcul, surtout dans les pre- 
miers instants de I'echange. 

Dans le cas ou Tarticle P est de forme complexe, on peut le decompo- 
ser en un ensemble de formes elementaires auxquelles on applique le maillage 
30 defini ci-dessus ou tout autre maillage adapte a la forme de I'article P. 
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Sur la figure 4, on a represents schematiquement I'enceinte de traite- 
ment 2 des articles alimentaires P telle que modelisee pour la mise en oeuvre de 
I'etape 22 de prediction de son comportement 

La prise en compte, dans le procede de determination des conditions 
de fonctionnement, des caracteristiques thermodynamiques et physiques de 
I'enceinte 2 lors de I'etape 22 de prediction du comportement de I'enceinte 2, est 
basee sur sa modelisation sous forme de tranches elementaires. 

Comme cela a ete decrit precedemment en reference a la figure 1, 
I'enceinte 2 de refroidissement est associee a un convoyeur 6. Elle est alimentee 
en fluide cryogenique 4 via a une ligne d'alimentation 5. L'enceinte 2 est assimi- 
lee a un parallelepipede rectangle. 

Pour determiner le profil theorique des temperatures du fluide 4, la me- 
thode mise en oeuvre pour la prediction du comportement de I'enceinte 2 consiste 
a effectuer une succession de bilans thermiques locaux. 

Dans ce but, on considere une modelisation du systeme thermodyna- 
mique de I'enceinte 2, en regime etabli, sous la forme de tranches elementaires 
34i a 34 n , perpendiculaires a sa longueur. La somme de ces tranches elementai- 
res 34^ a 34 n , represente le volume interne de I'enceinte 2. 

Pour chaque tranche elementaire 34i a 34 n , on effectue le bilan des 
transferts thermiques afin de determiner Tenthalpie du fluide 4 et done sa tempe- 
rature. 



Ce bilan doit tenir compte : 

- des deperditions thermiques avec Texterieur de I'enceinte 2 ; 

- du liquide cryogenique 4 injecte dans les zones de pulverisation ; et 

- des echanges entre les produits P et le fluide 4. 

Dans le cas de la tranche 34j de I'enceinte 2 de dimensions LTh, le bi- 



lan thermique se traduit par I'equation suivante : 





Dans cette equation : 

fj m correspond a Tenthalpie du fluide cryogenique 4 en sortie de la 
tranche elementaire 34j exprime en joules par kilogramme (J/Kg) ; 



WO 2004/092667 



PCT/FR2004/050121 



14 

H Mn corres P° ncl * I'enthalpie du fluide cryogenique 4 en entree de la 
tranche SISmentaire 34 s exprimee en joules par kilogramme (J/Kg) ; 

ff correspond a I'enthalpie liquide du fluide cryogenique 4 injects 

exprimee en joules. par kilogramme (J/Kg) ; 
5 Hmo corres P oncl a I'enthalpie de I'article P en entree de la tranche 

34j exprimee en joules par kilogramme (J/Kg) ; 

H M0 correspond a I'enthalpie de Particle P en sortie de la tranche 

34j exprimee en joules par kilogramme (J/Kg) ; 

]£ T correspond au coefficient d'echange thermique de I'enceinte 2 

10 avec I'extSrieur exprimS en watt par metre carrS et par kelvin (W/(m 2 K)) ; 

ihjsprayii) corres P oncl au debit massique de fluide cryogenique 4 vapori- 
se dans la tranche 34-, exprime en kilogramme par seconde (Kg/s) ; 

]jj Mi) correspond au debit massique de fluide cryogenique 4 entrant 

dans la tranche 34j exprime en kilogramme par seconde (Kg/s) ; 
15 jj2 p correspond au debit massique de produits a traiter exprime en ki- 

logramme par seconde (Kg/s) ; 

T Amb correspond a la temperature ambiante exprimee en kelvin ; et 

J 1 correspond a la temperature du fluide cryogenique 4 en entree de 

la tranche 34j exprimee en kelvin. 
20 Sur la figure 4 on a egalement represents les flux thermiques : 

ihjsr^HfL* est represents par la lettre A ; 

lh m H Mi) est represents par la lettre B ; 

ffl Mi) H m est represents par la lettre C ; 

lh p H mo est re P r ® sent6 P ar ia ,ettre D » et 
25 ihpHpsio est re P r esente par la lettre E ; 

aVeC Jfl M0 + A jh /Sprays 

Par experience, on sait que dans certaines conditions de fonctionne- 
ment (dSbit de production trop faible ou temperature du fluide cryogenique 4 trop 
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basse), le liquide cryogenique 4 injecte n'est que partiellement vaporise et une 
fraction du liquide s'ecoule vers I'entree de I'enceinte 2. 

Si Ton souhaite prendre en compte ce phenomene, il est preferable de 
resoudre les bilans locaux en commengant par la tranche elernentaire situee en 
5 sortie de I'enceinte 2. On fait done les calculs dans le sens inverse du trajet des 
articles P, soit selon I'axe X tel que represents sur la figure 4. 

Dans ce sens en effet, on peut reporter la fraction de liquide non vapo- 
rise dans la suivante et ainsi de suite jusqu'a atteindre les zones de ventilation ou 
les debits injectes sont nuls et ou les surplus liquides sont vaporises. 
10 Pour determiner la fraction de liquide cryogenique 4 non vaporise dans 

une tranche elernentaire, on designe une enthalpie du fluide limite, en-dessous 
de laquelle apparaTtra un titre liquide. 

Cela revient a fixer une temperature de fluide gazeux minimale dans 
I'enceinte 2. 

15 Le titre liquide non vaporise en sortie de la tranche elernentaire 34 j 

correspond a Xl(o et s ' ex P rim e sous la forme suivante : 

m Mo 

Si on simplifie les calculs en considerant que Tenthalpie de cette frac- 
tion liquide est sensiblement egale a I'enthalpie du fluide cryogenique 4 injecte, 
20 on obtient 1'expression du titre du liquide suivante : 

-H 

Dans cette equation, jum correspond a Tenthalpie limite de forma- 
tion d*un titre liquide dans une tranche elernentaire de I'enceinte 2. 

Sur la figure 5, on a represents le detail du cycle de test 29 et notam- 
25 ment de Tetape 20 de prediction de la temperature des articles P en sortie de 
I'enceinte 2. 

Afin de pouvoir realiser une prediction de la temperature des articles P 
en sortie de I'enceinte 2, le procede de refroidissement fait intervenir Petape 22 
de prediction du comportement de I'enceinte 2 et Tetape 24 de prediction du 
30 comportement des articles P. 
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On commence par realiser I'etape 22 de prediction du comportement 
de I'enceinte 2, lors d'une etape 40. 

Cette etape 40 delivre les pertes thermiques 42 par tranche elemen- 
taire, qui sont reintroduces dans I'etape 22. 
5 Apres avoir reitere cette operation un certain nombre de fois, on ob- 

tient les pertes thermiques totales 44 ainsi que le profil 46 des temperatures du 
fluide 4 dans I'enceinte 2, lors de I'etape 40. 

Pour calculer le bilan thermique de chaque tranche, I'etape 22 requiert 
les variations enthalpiques des articles P. De fait, lors de la premiere iteration, le 
10 profil des temperatures du fluide 4 dans I'enceinte 2 ne pouvant pas etre calcule, 
il est fixe arbitrairement. 

On met ensuite en oeuvre I'etape 24 de prediction du comportement 
des articles P, lors d'une etape 50. Cette etape 50 delivre I'enthalpie 52 de 
Tarticle P en sortie de I'enceinte 2, soit sa temperature. 
15 Eventuellement, I'etape 24 de prediction du comportement des arti- 

cles P delivre egalement les variations enthalpiques 54 d'un article P pour cha- 
que tranche elementaire de I'enceinte 2. Dans ce cas, cette information est re- 
tournee a I'etape 22 de prediction de comportement de I'enceinte 2, qui I'insere 
dans le bilan thermique de chaque tranche elementaire. 
20 L'enthalpie 52 de I'article P en sortie de I'enceinte 2 ainsi que le profil 

46 des temperatures du fluide 4 dans I'enceinte 2 et les pertes thermiques totales 
44, sont mis en relation afin de determiner le debit total du fluide, a I'etape 60. 

Eventuellement, on obtient egalement le debit 62 injecte dans chaque 
tranche elementaire. Dans ce cas, cette information est retournee a I'etape 22 de 
25 prediction du comportement de I'enceinte 2, qui I'insere dans le bilan thermique 
de chaque tranche elementaire. 

On verifie ensuite si le profil des temperatures du fluide 4 dans 
I'enceinte 2 est stable, a I'etape 70. 

Par exemple, le profil des temperatures du fluide est considere comme 
30 stable s'il repond deux fois de suite au critere suivant : 

H *^~ H ^ k -" <dif_profd 

H P s{\,k) 

Dans cette equation, difjprofil est une constante fixee par I'operateur. 
Au premier passage, le profil est considere comme instable. 
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Tant que le profil est considere comme instable, on retourne a I'etape 
40 et on recommence ta succession cooperations permettant de definir un profil. 

Une fois qu'un profil stable a ete obtenu, on verifie lors de I'etape 26, si 
la consigne determinee lors de I'etape 26 et portant sur la temperature des arti- 
5 cles P en sortie de I'enceinte 2 a ete atteinte. 

Si la consigne a ete atteinte, les parametres de fonctionnement abou- 
tissant au dernier profil des temperatures du fluide 4 a I'interieur de I'enceinte 2 
sont enregistres lors de I'etape 30 et forment une recette. 

Si la consigne n'a pas ete atteinte, a I'etape 28 les parametres de 
10 fonctionnement de Installation sont modifies. Par exemple, cette modification 
comporte une correction 102 sur le debit du fluide 4 avant de reiterer I'algorithme. 
Eventuellement, la modification comporte une correction 104 directement sur les 
parametres de fonctionnement conditionnant le profil des temperatures du fluide 
4, qui est utilisee a I'etape 22 de prediction du comportement de I'enceinte 2. 
15 Dans cet exemple, la prediction de la temperature des articles en sor- 

tie de I'enceinte est utilisee pour determiner le debit du fluide 4 injecte dans une 
enceinte cryogenique 2. 

On peut de la meme maniere influer sur le temps de sejour des articles 
P dans I'enceinte 2 en modifiant la vitesse du convoyeur 6 ou les temps d'arret 
20 dans le cas d'un convoyeur sequentiel. On peut egalement agir sur la vitesse 
d'extraction des gaz ou sur le taux de chargement, ou sur une combinaison de 
ces parametres. 

De meme, il est possible d'influer sur la balance entre les entrees d'air 
et les sorties de gaz, la vitesse d'extraction des gaz, la mise en vitesse des gaz, 
25 ou encore la recirculation des gaz en influant sur les elements de controle de ces 
parametres. 

Par ailleurs, si on met en oeuvre le procede de I'invention pour une ins- 
tallation disposant de capteurs sans contact de la temperature des articles en 
sortie, par exemple des capteurs bases sur le rayonnement thermique ou Timage 
30 infrarouge, ou encore sur une mesure par thermometrie micro-onde (TMO), tel 
que le capteur decrit dans le brevet FR-A-2 771 552, les resultats delivres lors de 
I'etape de prediction de la temperature des articles en sortie de I'enceinte peu- 
vent alors etre recoupes avec les mesures delivrees par ces capteurs. 
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Dans ce cas, on utilise Tune ou I'autre des informations pour verifier 

I'autre. 

Dans un autre mode de realisation, Tinformation delivree par le capteur 
est utilisee pour corriger la prediction. 
5 Bien qu'un mode particulier de realisation ait ete decrit, il n'est pas 

considere comme limitatif de la portee de la presente invention. 

En variante, le procede de refroidissement de Tinvention peut aussi 
etre applique dans une installation de froid mecanique presentant un dispositif 
d'echange de chaleur indirect. 
10 L'invention a ete decrite dans le cas du refroidissement d'articles ali- 

mentaires, cependant on peut Tappliquer egalement a d'autres types d'articles 
notamment metalliques. 

De plus, le terme refroidissement couvre egalement les systemes vi- 
sant au maintien et au controle d'une temperature inferieure a la temperature ini- 
15 tiale d'un article. 

Le procede de l'invention peut etre mis en oeuvre a I'aide par exemple 
d f un programme execute sur un ordinateur ou encore toute autre solution log i- 
cielle et/ou materielle adaptee. 
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REVEND1CATIONS 

1- Procede de determination des parametres de fonctionnement d'une 
installation de refroidissement thermique d'articles (P), comprenant une enceinte 
(2) traversee par lesdits articles (P) d'une entree a une sortie et utilisant un fluide 
5 de refroidissement (4), le procede comprenant : 

- une etape (16) de determination d'une consigne de temperature d'ar- 
ticles (P) en sortie de ladite enceinte (2) ; 

- une etape (18) de determination de parametres de fonctionnement 
initiaux pour ladite installation ; et 

10 - un cycle (29) de test des parametres de fonctionnement compre- 

nant : 

- une etape (20) de prediction de la temperature d'articles (P) en 
sortie de ladite enceinte (2) ; 

- une etape (26) de comparaison de la temperature de consigne 
15 avec la temperature predite ; et 

- si ladite etape (26) de comparaison revele un ecart superieur a 
un seuil predetermine, une etape (28) de modifications de parametres de fonc- 
tionnement de ladite installation (2) et une reiteration du cycle de test, 

ladite etape (20) de prediction etant realisee a partir de parametres de 
20 fonctionnement de ladite enceinte (2), de caracteristiques thermodynamiques et 
physiques de ladite enceinte (2) et de caracteristiques thermodynamiques et 
physiques desdits articles (P). 

2. Procede selon la revendication 1 , caracterise en ce que ladite etape 
(20) de prediction comporte une etape (22) de prediction du comportement de 

25 ladite enceinte (2) fondee sur la resolution de bilans thermiques sur des tranches 
elementaires du volume de ladite enceinte (2), realisee au moins a partir de ca- 
racteristiques thermodynamiques dudit fluide de refroidissement (4) et de carac- 
teristiques thermodynamiques et physiques de ladite enceinte (2). 

3. Procede selon la revendication 2, caracterise en ce que ladite etape 
30 (22) de prediction du comportement de ladite enceinte (2) est realisee en outre, a 

partir de parametres de fonctionnement de ladite installation. 

4. Procede selon la revendication 3, caracterise en ce que lesdits pa- 
rametres de fonctionnement de ladite installation represented au moins Tun des 
elements choisis dans le groupe consistant en : 
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- la vitesse (fun convoyeur (6) de transport desdits articles (P) au tra- 
vers de ladite enceinte (2) ; 

- le taux de chargement ; et 

- la ventilation de ['atmosphere de ladite enceinte (2). 

5 5. Procede selon Tune quelconque des revendications 1 a 4, characteri- 

se en ce que ladite etape (20) de prediction comporte une etape (24) de predic- 
tion du comportement desdits articles (P) fondee sur la resolution de ('equation 
de conservation de la chaleur discretisee et appliquee a un reseau de points spa- 
tiaux et temporels constituant un maillage desdits articles (P), realisee au moins a 
10 partir de caracteristiques thermodynamiques et physiques desdits articles (P). 

6. Procede selon la revendication 5, caracterise en ce que ladite etape 
(24) de prediction du comportement desdits articles (P) est realisee en outre, a 
partir de parametres de fonctionnement de ladite installation. 

7. Procede selon la revendication 6, caracterise en ce que lesdits pa- 
15 rametres de fonctionnement de ladite installation comprennent la temperature 

desdits articles (P) en entree de ladite enceinte (2). 

8. Procede selon Tune quelconque des revendications 5 a 7, caracteri- 
se en ce que ladite etape (24) de prediction du comportement desdits articles (P) 
est optimisee par des calculs de modification dudit maillage desdits articles (P) 

20 selon des suites mathematiques. 

9. Procede selon Tune quelconque des revendications 5 a 8, caracteri- 
se en ce que ladite etape (24) de prediction du comportement desdits articles (P) 
est optimisee par suppression des calculs de prediction pour des points spatiaux 
et temporels dudit maillage desdits articles (P) pour lesquels les variations 

25 d'enthalpie sont inferieures a un seuil predetermine. 

10. Procede selon les revendications 2 et 5 prises ensembles, caracte- 
risee en ce que ladite etape (20) de prediction de la temperature desdits articles 
(P) en sortie de ladite enceinte (2) est fondee sur ladite etape (22) de prediction 
du comportement de ladite enceinte (2) ainsi que sur ladite etape (24) de predic- 

30 tion du comportement desdits articles (P). 

11. Procede selon Tune quelconque des revendications 1 a 10, carac- 
terise en ce que ladite etape (28) de modification des parametres de fonctionne- 
ment comprend une etape de modification manuelle d'au moins une partie des 
parametres de fonctionnement. 



WO 2004/092667 



PCT/FR2004/050121 



21 

12. Procede selon Tune quelconque des revendications 1 a 11, carac- 
terise en ce que ladite etape (28) de modification des parametres de fonctionne- 
ment comprend la modification automatique d'au moins une partie desdits para- 
metres de fonctionnement. 
5 13. Procede selon Tune quelconque des revendications 1 a 12, carac- 

terise en ce que ladite etape (28) de modification des parametres de fonctionne- 
ment comprend la modification d'au moins un des parametres choisis par le 
groupe consistant en : 

- ie debit dudit fluide de refroidissement (4) ; 

10 - le temps de sejour desdits articles (P) dans ladite enceinte (2) ; 

- le debit de gaz extrait de ladite enceinte (2) ; 

- la mise en vitesse des gaz ; 

- la recirculation des gaz ; et 

- la balance entre les entrees d'air et les sorties de gaz. 
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